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Abstract 

We present predictions for the production cross section of a Standard Model 
Higgs boson in association with a tt pair at next-to-leading order accuracy using 
matrix elements obtained from the HELAC-Oneloop package. The NLO predic- 
tion was interfaced to the PYTHIA and HERWIG shower Monte Carlo programs with 
the help of POWHEG-Box, allowing for decays of massive particles, showering and 
hadronization, thus leading to final results at the hadron level. 

PACS 12.38.-t - Quantum ChromoDynamics 
PACS 13.87.-a - Jets in large-Q 2 scatterings 
PACS 14.65.Ha - Top quarks 
PACS 14.80.Bn - Standard Model Higgs boson 



1 Introduction 



In recent years a plethora of NLO QCD calculations have been presented in the literature. 
High-energy experiments, notably at the LHC, have and will be benefited by the progress 
in our computational ability to deal with higher order corrections to scattering amplitudes 
with many partons involved. In order though to get the optimum benefit and to produce 
predictions that can be directly compared to experimental data at the hadron level, a 
matching with parton showers and hadronization is ultimately inevitable. 

According to the literature, two methods exist that deal with the matching of parton 
showers (PS) to matrix element (ME) calculations at NLO accuracy, MC@NLO [Ul2j and 
POWHEG [3J H]. At the same time, several computational frameworks for NLO QCD 
calculations have emerged. In constructing a general interface of PS to ME computations 
with NLO accuracy, we have chosen to combine the POWHEG and HELAC-NLO [5] 
approach, respectively. The first application of our project was presented in Ref. [6]. 



In this letter we continue our investigation and validation of such an interface by 
studying Higgs boson production in association with a tt pair. There are several reasons 
to perform this work. On the phenomenological side, there is still interest to use this 
channel to search for Higgs boson in the low mass regime [7]. On the theoretical side, 
the appearance of recently developed tools as implemented in aMC@NL0 [8], also used 
in the analysis of Higgs boson production in association with a tt pair, offers a unique 
opportunity for cross-checking and validating our computations and the corresponding 
computer programs. 



We performed our calculations using the POWHEG-Box [9j, a flexible computer framework 
implementing the POWHEG method, already used in the study of different processes also 
by other authors (among the others, the most recent applications include Ref. [TOj [TIJ [12J 
[13]). In this work the following ingredients were requested and provided to POWHEG-Box 
as input: 

• The flavor structures of the Born (gg— > Htt, qq— > Htt, qq— > Htt) and real radiation 
emission (qq— > Httg, gg—¥ Httg, qg— > Httq, gq— > Httq, qq— > Httg, qg— > Httq, gq— > 
Httq) subprocesses (qG{u,d,c,s,b}) were specified. 

• The Born-level phase space was obtained by using the invariant mass of the tt quark 
pair and four angles. 

• Squared ME's with all incoming momenta for the Born and the real-emission pro- 
cesses were built using amplitudes computed by HELAC-Oneloop [TJ], on the basis 
of skeleton structures generated by this same program and HELAC-PHEGAS [15] , re- 
spectively. The ME's in the physical channels were obtained by crossing. 

• Color structures for building the color-correlated Born amplitudes were taken from 
HELAC-Dipoles [16] . 

• Spin-correlated Born amplitudes were projected from the helicity basis to the Lorentz 
one by using the polarization vectors. 
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• HELAC-Oneloop was also used for computing the finite part of the virtual correction 
contributions in dimensional regularization, on the basis of the OPP method [T7] 
complemented by Feynman Rules for the computation of the QCD R2 Rational 
terms [18]. 

With this input POWHEG-Box was used to generate events at the Born level plus first 
radiation, stored in Les Houches Event Files (LHEF) [19]. Then, one can choose any 
shower Monte Carlo (SMC) program for generating events with hadrons, just by showering 
the previous ones and applying hadronization and decay models (it is worth mentioning 
that POWHEG-Box was written in the effort of being shower detail independent). We 
worked with both the PYTHIA [201 El] and HERWIG [221 E3] SMC generators. 

In the POWHEG-Box the first emission is the hardest one measured by transverse mo- 
mentum which is also the ordering variable in PYTHIA. However, if the ordering variable in 
the shower is different from this one, such as in HERWIG, the hardest emission can be dif- 
ferent from the first one. In such cases, in order to face this issue, HERWIG is constrained 
to discard shower evolutions (vetoed shower) with larger transverse momentum in all 
splittings occuring after the first emission, that is computed by POWHEG-Box . In addi- 
tion, a truncated shower, simulating wide-angle soft emission before the hardest emission, 
is in principle also needed for POWHEG matching [3], but its effect turned out to be 
small [21]. As there is no implementation of truncated shower in HERWIG using external 
LHEF, the effect of the truncated shower is absent from our predictions. 

3 Checks 

The Born and the real ME's were tested in randomly chosen phase-space points against re- 
sults from HELAC-PHEGAS . The consistency among real emission, Born, color-correlated 
and spin-correlated ME's was tested in randomly chosen real-emission phase-space re- 
gions, by examining the ratio between the real-emission amplitudes in soft and collinear 
limits and the corresponding subtraction terms, automatically computed by POWHEG-Box 
as approaching the related singularities. 

The process presented here was studied extensively at the NLO accuracy in the liter- 
ature [23 ES, [23 EHl EH], which enabled us to make detailed checks of our calculation. 
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Table 1: Comparison of the NLO cross sections by POWHEG-HELAC (PH) to the predictions 
of Ref. [26] (lit) at the default scale (jlr = /xp = /jlq. 
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Figure 1: Transverse momentum and rapidity distributions of the Higgs boson at the 
Tevatron from our calculation by means of P0WHEG-HELAC (POWHEG) and from Ref. [26] 
(pure NLO). The lower panels show the ratio of the results (POWHEG/NLO). 



As for the virtual contribution, we compared our amplitudes to the predictions of 
MADL00P using the input parameters and phase space point specified in Ref. [SU] for 
different renormalization scale choices [31J, and in all cases we found agreement up to 5-6 
digits. 
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We then compared the cross-sections at the NLO accuracy to the predictions of 
Ref. [26]. At this purpose, we used the MRST2002nlo [32] PDF set, provided by LHAPDF 
[55] . with a 2-loop running a s , 5 light flavours and Af s = 239 MeV, and the following 
parameters: m# = 120 GeV, m t = 174 GeV, Gf = 1-16639 • lCT 5 GeV~ 2 . The renormal- 
ization and factorization scales were chosen equal to /i, where /i was varied in the /io/4 - 
4/xo interval, with default scale being fio = m t + ^i?/2. The calculations were performed 
at the planned LHC energy of \fs = 14TeV. The results of this comparison are presented 
in Table [TJ showing agreement. We also made a comparison to the predictions of Ref. [29] , 
using the CTEQ5M PDF set from LHAPDF, a 2-loop running a s with a s (M z ) = 0.118, m t 
and Gf as given above, and again found agreement within their quoted uncertainty. 

Finally, we also compared differential distributions from events already including first 
radiation emission with the NLO predictions of Ref [26J. We found agreement for all 
distributions presented in that paper. As examples, we show in Fig.[T] the transverse 
momentum and rapidity distributions of the Higgs boson in the Tevatron setup R 



4 The effect of the shower 

In interfacing NLO calculations to SMC programs we mainly aim at estimating the effects 
of showers and hadronization, therefore, we analyzed the events at two different stages 
of evolution: 

Decay: starting from the events collected in LHEF produced by POWHEG-HELAC, 
we just included on-shell decays of t-quarks and the Higgs boson, as implemented 
in PYTHIA, and further decays of their decay products, like charged leptons (the r 
is assumed to be unstable), W and Z, turning off any shower and hadronization 
effect. 

Full SMC: decays, showering evolution and hadronization have been included in 
our simulations, using both PYTHIA and HERWIG. 

In both SMC setup muons (default in PYTHIA) and neutral pions were assumed as stable 
particles. All other particles and hadrons were allowed to be stable or to decay according 
to the default implementation of each SMC. Quark and Higgs masses, as well as W, Z 
masses and total decay widths, were tuned to the same values in PYTHIA and HERWIG. 
On the other hand, each of the two codes was allowed to compute autonomously partial 
branching fractions in different decay channels for all unstable particles and hadrons. 
Multiparticle interaction effects were neglected (default in HERWIG ). Additionally, the 
intrinsic spreading of valence partons in incoming hadrons in HERWIG was assumed to 
be 2.5 GeV. Considering this setup, we found agreement between PYTHIA and HERWIG 
predictions within 5 %, except in bins where the statistics is very small. Beside the 
conceptual differences in the parton shower and hadronization algorithms between the 
two SMC generators, written on the basis of different theoretical ideas (p± vs. angular 
ordering, string model vs. cluster hadronization and preconfinement), a possible origin 
of this overall small discrepancy is the absence of the truncated shower in the HERWIG 

1 Note that these comparisons refer to distributions without normalization, as opposed to the plots 
in the original publication. 
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prediction. Unfortunately, however, as already mentioned, we cannot check the last point 
within the POWHEG-Box framework. In any case, the modest entity of the discrepancy 
means that the effect of the truncated shower, not included in our analysis, is small. 
Further explanations of small discrepancies can also be related to the different decay 
channel scheme and branching fractions in the two SMC codes. 

In our computation, we adopted the following parameters: \/s = 7 TeV, CTEQ6 . 6M 
PDF set from LHAPDF, with a 2-loop running a s , 5 light flavours and A^ s = 226 Me V, 
m t = 172 GeV, m H = 120 GeV , G F = 1.16639 • 10" 5 GeV~ 2 . The renormalization and 
factorization scales were chosen equal to /j,q = m t + m#/2. We decided to switch on all 
possible decay channels of the Higgs boson, implemented in the SMC programs]^] We 
used the last version of the SMC codes: PYTHIA 6.425 and HERWIG 6.520. 

We studied the effect of the full SMC by comparing distributions at the decay and 
SMC level. As the number of particles is very different at the end of the two stages, we 
first made such a comparison without any selection cut, in order to avoid the introduction 
of any bias. This way the cross-section at all levels is indeed exactly the same. We found 
Cpowheg-helac = °"powheg-helac+decay = °"pqwheg-helac+smc = 95.872 ± 0.007 fb (here and in 
all following a predictions the quoted uncertainties are the statistical ones only). As 
an illustrative example, we present the distributions of the transverse momentum and 
rapidity of the hardest jet, p±j 1 and y^, in Fig. [2] and Fig.|3j respectively. The jets are 
reconstructed through the anti — k± algorithm with R = 0.5, by using Fast Jet 2.4.3 [34J. 
One can observe a rather significant softening in the transverse momentum spectrum as 
going from results at the decay level to full SMC ones. On the other hand, the effect of 
the shower on the rapidity of the hardest jet is almost negligible and rather homogeneous. 
We distinguished several classes of jets including final state emissions, according to their 
origin: jets that can be traced back to (i) first radiation emissions, (ii) the decay products 
of the Higgs boson, (iii) the decay products of the top and antitop quarks, (iv) a mixing 
of the previous ones. In particular, main contributions to the Pi_,ji spectrum shown in 
Fig. [2] are due to jets of the (iv) and (iii) class. As for the y^ distribution, the tails are 
dominated by jets of the (iv), (i) and (iii) classes. 

In Ref. [30] the invariant mass, vtibbi and the separation in the rapidity-azimuthal- 
angle plane of the two hardest lowest-lying B-hadrons, AR BB , were studied, by choosing 
a dynamical scale for the generation of the hard-scattering events and by taking into 
account only the H — > bb decay channel. Besides computing all pt distributions already 
presented in that paper, always finding agreement, we also computed the aforementioned 
B-hadron distributions reproducing the same simulation setup and without applying any 
cut, considering only the H — > bb decay channel, as well as all channels. In the former, we 
found agreement with the predictions of Ref. [30J . The effect of the remaining channels, 
not studied in that work, produces an increase in the region below 80 GeV in the ttibb 
spectrum and only for large ARbb, as shown in Figs. [4]and[5j Pairs of B-hadrons both 
including quarks that can be traced back to Higgs decays only populate the region of the 
rriBB spectrum below uih- The region above m# is instead dominated by pairs with at 
least one B-hadron that can be traced back to a (anti)top decay. 

The effects of the decay and shower also depend on the selection cuts. While the 
typical selection cuts include both leptons and hadrons, we started with a restricted 

2 In PYTHIA there are two more decay channels than in HERWIG, and partial decay fractions in each 
leptonic, bosonic and partonic channel differ in the two codes. 



5 



> 

CD 

o 



10" 



10 - 



10" - 



3 



10" 



10 

o 1.5 E" 
1.0 E 



I 1 I 1 I 1 

Decay (PH+PYTHIA) 
SMC (PH+PYTHIA) 



=■=4 




-ln ff = 120GeV 
fj, = m t + m H /2 
CTEQ6 . 6M 

anti-fcx, ii = 0.5 

III I 



,_ r • 



0.5 E- |- 
O.o F-T I 



1 



1 



1 



1 



J 3 



100 200 300 400 

P±,h [GeV] 



500 600 



Figure 2: Transverse momentum distribution of the hardest jet. The results at the decay 
level (dashed line) are compared to the ones at the SMC level (solid line), obtained by 
interfacing P0WHEG-HELAC to PYTHIA . No selection cuts were applied. 
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Figure 3: The same as Fig. [2j as for the rapidity of the hardest jet. 



6 



10 m | i | i | i | i | i | i | i | i | ig 



2Q 



2 

io" 5 

o 1.5 
1.0 
0.5 
0.0 



= 7TeV 
rn t = 172.5GeV 
MSTW2008: 



/' 



(m± t ■ m I t ■ m_L, H ) 1/3 




H-» bb channel (PH+HERWIG) 

all channels (PH+PYTHIA) 

all channels (PH+HERWIG) 

I I I I I I I I I I I I I I I I I J 



'MMlVir' 



P-u-j— ■ 

' g 



20 40 60 80 100 120 140 160 180 200 

m B B [GeV] 



Figure 4: Invariant mass distribution of the two hardest lowest-lying B hadrons at the 
SMC level. These results are presented in sigma per bin, just to allow an easier comparison 
with the results of Ref. [30], having been obtained in the same setup, without applying 
any cut. The effects of including all H decay channels, with respect to the case of a 
single H — > bb channel (dashed line), were computed by interfacing P0WHEG-HELAC to 
both PYTHIA (solid line) and HERWIG (dash-dotted line). 



set, not involving any leptonic cut, in particular as for the number of leptons. This is 
motivated by the fact that this number can be quite different at the decay and shower 
level, as the shower produce many secondary leptons. Thus, coming back to our setup, 
we considered and implemented cuts only on jet variables: (i) pj_ min = 20 GeV and (ii) 
< 2.5 for all jets, (iii) #jets > 4 in each event. We show the distribution of the 
scalar sum of all transverse momenta in the event, H± in Fig. [6} We can see that the 
spectrum becomes softer due to showering effects, with respect to the one computed at 
the decay level, as can be understood on the basis of PS physics. The same is true and 
even more evident if one singles out the hadronic component of H± , as well (not shown) . 
On the other hand, the effects of the shower on the (anti-)lepton transverse momentum, 
pj_, and the missing transverse momentum, p™ lss , as shown in Figs. [7] and |8j are small 
and rather uniform, except for a significant increase for small values, which is due to 
secondary leptons produced in the shower. The cross-sections at the decay and at the 
SMC level, after the hadronic cuts listed above, amount to o>owheg-helac+decay = 92.29 ± 

0.01 ft), (TpQWHEG-HELAC+PYTHIA = 90.46 ± 0.01 ft and OPOWHEG-HELAC+HERWIG = 90.99 ± 0.01 ft. 
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Figure 5: The same as Fig. |4| as for the azimuth-rapidity distance correlations of the two 
hardest lowest-lying B hadrons. 



5 Predictions 

In this section we present predictions for Htt production with parton shower and hadroniza- 
tion effects at LHC. In addition to the jet cuts (i-iii) mentioned above, we also applied 
selection cuts on the leptonic variables]^] (iv) we focused on the dileptonic channel, asking 
for exactly one £ + and one t~ with (v) p e ± min = 20 GeV and (vi) \y e \ < 2.5, whereas the 
transverse missing energy of the event was constrained to (vii) -$_i_ jmin > 30 GeV. 

In Fig.[9]we present the distribution of the invariant mass of all jet pairs. Here the 
effect of the shower is again quite significant. In particular, there is a small bump around 
the Higgs mass, as already noticed in the literature (35], visible in the data at the decay 
level, which is completely washed out when PS is included. The same is true for the 
invariant mass distribution of the two hardest jets (not shown). 

The cross-section after cuts, at the SMC level, were found to be o" P o WH eg-helac+pythia = 
5.376 ± 0.010 fb and ctpowheg-helac+herwig = 5.521 ± 0.011 fb. 

One of the biggest differences in the results produced by PYTHIA and HERWIG inter- 
faced to the POWHEG-Box noticed in Ref. [12] in the study of a different process was the 
observation that HERWIG gives rise to hard jets more central than PYTHIA. We observe the 
same trend in our results, but by far to a lesser extent. In particular, in the bins around 
zero rapidity the ratio between the rapidity distributions of POWHEG-HELAC + HERWIG and 
POWHEG-HELAC + PYTHIA found in our study amounts to maximum 1.05, both in case of 
the hardest and the second hardest jet. On the other hand, the agreement between the 
two SMC, as for the rapidity distributions of leptons and antileptons, was found to be 

3 Similar cuts are applied by the LHC experiments. 
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Figure 6: Distribution of the sum of the transverse momenta of all particles, at the decay 
level (dashed line) and at the SMC level, as obtained by interfacing P0WHEG-HELAC to 
PYTHIA (solid line) and HERWIG (dash-dotted line). Only hadronic cuts were applied. 
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Figure 7: The same as Fig. [6j as for the transverse momentum of all antileptons. 
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Figure 8: The same as Fig. [6j as for the total missing transverse momentum. 
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Figure 9: Invariant mass distribution of all jet pairs. Results at the decay level (dashed 
line) are shown together with results after showering and hadronization, obtained by 
interfacing P0WHEG-HELAC predictions to both the PYTHIA (solid line) and HERWIG (dash- 
dotted line) SMC. The full set of hadronic and leptonic cuts was applied. 
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6 Conclusions 

We would like to conclude this letter by emphasizing that interfacing NLO calculations, 
as structured in HELAC-NLO, with PS and hadronization effects, within the POWHEG 
framework, open the exciting possibility of realistic, precise and reliable simulations 
of pp hard-scattering collisions, leading to multiparticle production, up to the hadron 
level, as detected in the accelerator experiments. We have presented predictions for 
the signal, pp — > ttH production, illustrating their consistency with existing calcula- 
tions, as well as the phenomenological potential of our approach. A thorough analy- 
sis, including the most important backgrounds, notably pp — > ttbb, within the same 
framework, is in progress and will be reported elsewhere. The event files produced by 
POWHEG-HELAC, together with further detail and results of our project, are available at 
http : //grid. kf ki .hu/twiki/bin/view/DbTheory/TthProd. 
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